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1^ ■ Based on a sample corresponding to 4.3 x 10^ produced r-pair events, we have 

'^ ■ studied hadronic dynamics in the decay r~ — > Ut-tt^'it^tt^ in data recorded 

by the CLEO II detector operating at the CESR e+e" collider. The decay 
^~~^ [ is dominated by the process r~ — > ly-rCii (1260) , with the a^ meson decay- 

pg I ing to three pions predominantly via the lowest dimensional (mainly 5-wave) 

CN ' ai -^ p^TT^ Born amplitude. From model-dependent fits to the Dalitz plot 

~j ■ and angular observables in bins of Svr mass, we find significant additional con- 

O ■ tributions from amplitudes for oi decay to air, /o(1370)7r and /2(1270)7r, as 

^ i well as higher dimensional oi — > /ovr and p'lr amplitudes. Notably, the squared 

■"^ ' an amplitude accounts for approximately 15% of the total t~ — > UrTT'TT^ir^ 

Q ' rate in the models considered. The data are well described using couplings 

Q^' to these amplitudes that are independent of the Svr mass. These amplitudes 

D . also provide a good description for the r~ -^ I'T-n" tt^ tt" Dalitz plot distribu- 

. . i tions. We have searched for additional contributions from t~ -^ ^'T-7r'~(1300). 

. !^ I We place 90% confidence level upper limits on the branching fraction for this 

p\ ' channel of between 1.0 x 10^^ and 1.9 x 10~^, depending on the vr' decay mode 

j^ ' considered. The vr^vr'^vr'^ mass spectrum is parametrized by a Breit-Wigner 

form with a mass-dependent width which is specified according to the results 
of the Dalitz plot fits plus an unknown coupling to an oi -^ K*K amplitude. 
Prom a x^ ht using this parametrization, we extract the pole mass and width 
of the ai, as well as the magnitude of the K*K coupling. We have also in- 
vestigated the impact of a possible contribution from the a'^(1700) meson on 
this spectrum. Finally, exploiting the parity-violating angular asymmetry in 
ai — > Svr decay, we determine the signed value of the r neutrino helicity to be 
h^^ = — 1.02ib0.13 (stat.)ib0.03 (syst.+model), confirming the left-handedness 
of the r neutrino. 
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I. INTRODUCTION 

The decay t~ ^ i/T-[37r]~ has been the subject of much interest over the years. Because of 
the transformation properties of the weak current under parity and G-parity, r lepton decay to 
an odd number of pions is expected to occur exclusively through the axial vector current, ignoring 
isospin-violating effects. Thus the Svr system in this decay must have spin-parity quantum numbers 
J = 0~ or 1"*". As a result of this plus the purely weak interaction involved in r decay, such 
decays provide an excellent opportunity to investigate the axial vector hadronic weak current and 
the dynamics of axial vector meson decay. 

The T~ — > ^'T-[37r]^ decay is dominated by production of the poorly understood ai(1260) meson, 
which is believed to decay mainly via the lowest dimensional Born (mostly 5- wave) pvr intermediate 
state. The world average values ffl for its mass and width are 1230 it 40 and 250-600 MeV respec- 
tively, determined primarily from r decay. The theoretical understanding of the ai is not rigorous 
- many models have been proposed [0-0] to describe the line shape and resonant substructure, but 
none have provided an entirely satisfactory description of the data. Additional experimental input 
is essential for a better understanding of this system. 

Recent experimental studies of the r~ -^ z^^vr^vr+Tr^ decay have been carried out by the AR- 
GUS [|J|, OPAL @ and DELPHI |ll|,|l2| collaborations. Analyses of earlier data such as that 



by Isgur, Morningstar and Reader M had demonstrated the presence of D-w&ve pir production. 
With ~ 7, 500 events, ARGUS H measured the ratio of amplitudes at the nominal oi mass to be 



D/S = -0.11 ± 0.02. In a sample of ~ 6, 300 events, OPAL [|T0| found D/S = -0.10 ± 0.02 ± 0.02. 
Another ARGUS analysis Q of ~ 3, 300 lepton-tagged events considered many additional ampli- 
tudes. This analysis found a signal at the 4.2 standard deviation (a) level for the presence of 
an ai -^ /2(1270)7r amplitude. Neither ARGUS nor OPAL found evidence of non-axial- vector 
contributions such as production of the J^ = 0^ 7r'(1300) which decays via pir and crvr. 

In addition, the decays r^ -^ ^'T-[37r]^ and r^ -^ i/t-[5tt]~ have been employed to constrain the 



r neutrino mass |11,1^, through investigation of the endpoint in the invariant mass and energy 



spectra of the multi-pion system. Notably, ALEPH |14| has obtained an upper limit of 18.2 MeV 



at the 95% confidence level (CL) on the Uj- mass, based on these decays. These analyses rely on 
an understanding of the hadronic dynamics. DELPHI [OJl^, with ~ 6,500 events, has reported 
anomalous structure in t~ -^ Ut-tt^ tt^ n^ in the context of a Vr mass analysis. In that work, the 
Dalitz plot distribution for events with very high 3-7r mass (> 1.5 GeV) is suggestive of enhanced 
Z)-wave yOTT production, while the 2>ti mass spectrum shows an excess in this region relative to 
expectations from a single oi resonance. Inclusion of a radial excitation (an a'l) with a large D- 
wave coupling to pvr provides an improved description of the DELPHI data, but weakens the Vr 
mass limit. 

In this article, we present results from a model-dependent analysis of the decay r~ — > iyrT^~T:^T:^ , 
based on data collected with the CLEO II detector. We perform fits to models to characterize both 
the substructure as seen in Dalitz plot and angular variables, as well as the ai resonance parameters 
as seen in the 37r invariant mass spectrum. The 7r~7r^7r'' channel has several advantages relative to 
the all-charged mode. First, the multihadronic and r feed-across backgrounds are smaller. Second, 



^Generalization to charge conjugate reactions and states is implied throughout, except as noted. 



the 7r~7r''7r'' mode may be better suited for discerning substructure involving isoscalar mesons 
because there is only one pairing of pions which can have / = 0, unlike in the all-charged case. 
This second point is particularly relevant in light of the ARGUS result for f2ir production M and 



the recent observation of the decay r~ -^ z^T-/i(1285)7r~ by CLEO ||l^. These results suggest that 
isoscalars may play a role in ai decay. 

The outline for the article is as follows. In Section ||, we describe the data sample and event 
selection. We describe the basic elements of the model used to characterize the r~ -^ i't-it~tt^it^ 



data in Section III. In Sections |^, ^ and VI, we describe the three analyses of the hadronic 



structure: (1) performing fits to the substructure based on Dalitz plot and angular observables; (2) 
extending these fits to determine the signed r neutrino helicity; and (3) performing fits to determine 
the resonant structure of the Svr mass spectrum, making use of results from the substructure fits. 



We summarize the results and conclude in Section VII. 



II. DATA SAMPLE AND EVENT SELECTION 

The analysis described here is based on 4.67 fb^^ of e^e~ collision data collected at center-of- 
mass energies 2£'beam of ~ 10.6 GeV, corresponding to 4.3 x 10^ reactions of the type e~^e~ —>■ t^t^ . 
These data were recorded at the Cornell Electron Storage Ring (CESR) with the CLEO II detec- 
tor 116] between 1990 and 1995. Charged particle tracking in CLEO II consists of a cylindrical 



six-layer straw tube array surrounding a beam pipe of radius 3.2 cm that encloses the e~^e~ inter- 
action point (IP), followed by two co-axial cylindrical drift chambers of 10 and 51 sense wire layers 
respectively. Barrel (| cos9\ < 0.81, where 9 is the polar angle relative to the beam axis) and end 
cap scintillation counters used for triggering and time-of-flight measurements surround the tracking 
chambers. For electromagnetic calorimetry, 7800 CsI(Tl) crystals are arrayed in projective (toward 
the IP) and axial geometries in barrel and end cap sections respectively. The barrel crystals present 
16 radiation lengths to photons originating from the IP. 

Identification of r~ -^ Vt-tt'tt^t^^ decays relies heavily on the segmentation and energy resolution 
of the calorimeter for reconstruction of the vr'^'s. The central portion of the barrel calorimeter 
(I cos 9\ < 0.71) achieves energy and angular resolutions of as/E (%) = 0.35 E^-''^ + 1.9 — 0.1 E and 
o"0 (mrad) = 2.8 /\/E + 2.5, with E in GeV, for electromagnetic showers. The angular resolution 
ensures that the two clusters of energy deposited by the photons from a vr'^ decay are resolved over 
the range of vr'^ energies typical of the r decay mode studied here. 

The detector elements described above are immersed in a 1.5 Tesla magnetic field provided by 
a superconducting solenoid surrounding the calorimeter. Muon identification is accomplished with 
proportional tubes embedded in the flux return steel at depths representing 3, 5 and 7 interaction 
lengths of total material penetration at normal incidence. 

A. Event selection 

To identify events as tt candidates we require the decay of the r"*" (denoted as the 'tagging' 
decay) that is recoiling against our signal r~ decay to be classified as VrC'^Ue, Fr/W^f^, T^t'^'^, or 
FT-7r"''7r'^. Thus, we select events containing two oppositely charged barrel tracks separated in angle 
by at least 90°. To reject backgrounds from Bhabha scattering and two-photon interactions we 
require track momenta to be between 0.08£'beam and 0.90-E'beam- Clusters of energy deposition 



in the central region of the calorimeter (|cos0| < 0.71) that are not matched with a charged 
track projection are paired to form vr*^ candidates. These showers must have energies greater than 
50 MeV, and the invariant mass of the photon-pair must lie within 7.5a of the vr'^ mass where a 
varies between ~ 4 — 7 MeV. Those tt" candidates with energy above 0.06 -Bbeam after application 
of a vr*^ mass constraint are associated with any track within 90°. 

A vr^vr'^Tr'^ candidate is formed from a track which has two associated vr'^ candidates as defined 
above. If more than one combination of n^ candidates can be assigned to a given track, only one 
combination is chosen: namely, that for which the largest unused barrel photon-like cluster in the 
vr^vr'^Tr'' hemisphere has the least energy. A cluster is defined to be photon-like if it satisfies a 1% 
confidence level cut on the transverse shower profile and lies at least 30 cm away from the nearest 
track projection. 

As mentioned earlier, to reject background from multihadronic {e~^e~ — > qq) events, the tag 
system recoiling against the it~tt^it^ candidate must be consistent with r decay to neutrino(s) plus 
e"*", //"*", Tr"*" or tt~^it^ (denoted as p~^). The recoiling track is identified as an electron if its calorimeter 
energy to track momentum ratio satisfies 0.85 < E/p < 1.1 and if its specific ionization in the main 
drift chamber is not less than 2a below the value expected for electrons. It is classified as a muon if 
the track has penetrated to at least the innermost layer of muon chambers at 3 interaction lengths. 
If not identified as an e or a /x, then if the track is accompanied by a third n^ of energy > 350 MeV 
which lies closer to it than to the irlir^ system, the track-vr*^ combination is classified as a p tag. 
The invariant mass of this system must be between 0.55 and 1.20 GeV. If not identified as an e, 
/i, or p tag, the recoil track is identified as a single vr tag. To ensure that these classifications 
are consistent with expectations from r decay, events are vetoed if any unused photon-like cluster 
with I cos^l < 0.95 has energy greater than 200 MeV, or if any unmatched non-photon-like cluster 
has energy above 500 MeV. The missing momentum as determined using the ttItt^ and tagging 
systems must point into a high-acceptance region of the detector (| cos^missl < 0.9), and must have 
a component transverse to the beam of at least 0.06£'beam- 

B. Final event sample 

After all cuts, the remaining sample consists of 51,136 events. The normalized invariant masses 
of the two photon-pairs, Syy = {Myy — m^o)/a^^, are plotted against each other in Fig. |l](a). In 
Fig. |l|(b), S^^ is plotted for all n^ candidates along with the corresponding distribution from the 
T Monte Carlo (MC) sample described in the following section. We define the ir^n^ signal region 
to be that where —3.0 < S^^ < 2.0 for both n^ candidates. In the vr^vr*^ signal region there are 
36710 events, of which 17234 are tagged by leptonic decays of the recoiling r. To estimate the 
contributions from fake tt'^'s, we also define side and corner band regions using —7.5 < S^^ < —5.0 
and 3.0 < S.yy < 5.5. 

In Fig. |2|, we plot the 7r~7r'^7r'^ mass for events in the tt^tt'^ signal and side band regions, for data 
and r MC samples. The events above the r mass in Fig. |2|(c) are dominantly due to feedacross 
from T ^ pv decays where the second vr'^ is being picked up from the recoil r decay. The r Monte 
Carlo simulation accounts for most of this high mass tail, but not all, with the remainder being 
due to a small qq background contribution. The high mass events are essentially absent from the 
lepton-tagged events, plotted in Fig. ||(d). 
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FIG. 1. (a) Plot of the normalized photon-pair invariant mass S. 



m^o)/a^.y for the 



two vr'' candidates after all other cuts have been applied, (b) Comparison of the S^^ distribution 
(two entries per event) for data (points) and tt Monte Carlo (histogram) samples. 
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FIG. 2. ir^TT^iT^ mass spectrum from the all-tag (a) and lepton tag (b) samples, after cuts. 
Shown are events in the vr'^vr'^ signal (dark points for the data, unshaded histogram for the Monte 
Carlo spectrum) and side band (light points for the data, shaded histogram for the MC spectrum) 
regions. Plots (c) and (d) show the same spectra as in (a) and (b) plotted over a larger range and 
on a logarithmic scale. 



C. Monte Carlo samples 

For determination of detection efficiency for our signal decay as well as for backgrounds from 
other T decay modes, we rely on a sample of Monte Carlo e~^e~ -^ t~^t~ events with equivalent 
luminosity approximately three times that of the data. These events were generated using the 
KORALB/TAUOLA program, and then passed through the GEANT-based [jl8| CLEO II detector 



simulation package. The full CLEO event reconstruction program was then run on this sample. 
The MC distributions shown in Figs, ffl and |2| are derived from this sample. 

In TAUOLA, the r^ -^ ^'T-[37r]^ decay is described with a single ai resonance decaying solely via 
the lowest dimensional (s-wave, in the notation introduced in the next section) Born amplitude 
for (p + p')7r production, following the model of Kiihn and Santamaria [^. We have tuned the 
oi mass and width to yield a Stt mass spectrum that roughly matches that seen in our data in 
the all-charged mode. Although the data and MC Svr mass spectra show reasonable agreement on 
average (see Fig. |2|), close inspection reveals significant deviations, particularly in the high mass 
region located roughly between 1.4 GeV and the r mass. The Dalitz plot distributions agree poorly 
with the corresponding MC distributions, especially in the M^o^o projection and most strikingly 
at high Stt mass. 



For the substructure fits described in Section IV, we generated additional MC samples for 
our signal mode plus key background r decay modes. For these samples, we developed special 
purpose event generators. Unlike the treatment in KORALB/TAUOLA, we implemented radiative 
effects according to an approximation in which they factorize from the rest of the differential 



matrix element, as required by the reverse Monte Carlo approach described in Section IV A, In 
addition, the signal mode was generated with a Svr mass spectrum weighted towards high values 
so as to ensure high statistics in the high-mass region where the data show the most apparent 
deviation from the model used by TAUOLA in both Dalitz plot and Svr mass distributions. 

III. MODEL OF T- -^ i/^7r-7r°7r° 

Tau lepton decay to neutrino plus three pions follows the form 

^$4, (1) 



dr(r -^ UrSir) 



2mr 






where L^"^ represents the lepton tensor for weak decay, J'^ denotes the hadronic weak current for 
production of three pions, and d^4 is the Lorentz-invariant 4-body phase space element for the 
decay. The goal of this analysis is to probe the structure of the hadronic current, benefitting from 
the well-understood properties of the weak interaction. 

In principle J^ is comprised of vector and axial vector currents: 

,P = J^ + J'X =< 3^|d7^M|0 > + < 37r|d757;,u|0 >, (2) 

however G-parity conservation requires that J^ = 0. Thus we consider only contributions from the 
axial vector current. 

In r decay, the squared momentum transfer s = M^^ is small, and thus the dynamics are 
expected to be dominated by resonance production. The possible axial vector (J^ = 1+) resonance 
contributions are the ai(1260) and radial excitations, i.e., the a[. In addition, pseudoscalar (J^ = 

9 



0~) contributions are possible, i.e., the 7r'(1300), although these are expected to be suppressed 
according to the Partially Conserved Axial Current (PCAC) hypothesis. In this section we describe 
the model used to parametrize J^, assuming the Sir system is in a J^ = l"*" state. 

A. Model for substructure in 1+ — y Stt 



The strong decay of the ai is expected to involve substructure which is again dominated by 
resonance production. We write for the contribution to J'^ involving ai production 



Ja', 



Ba,{s)xJ2/3dt, 



(3) 



where i3a^(s) denotes the ai Breit-Wigner, /3j are complex coupling constants, and j^ contain form 
factors describing components of the substructure involving specific resonances. The details of this 
parametrization are given in Appendix ^. As an example, in the case of s-wave pir production. 



■/I 

Ipn, s-wave 



T^^'^ [qiuB,isi)FR^^{h)-q2uBpis2)FR^^ik2) 



(4) 



We define pi, p2 and p3 as the four-momenta of the three pions, in our case ttj*, vTg and vr", 
respectively. We define a = pi + p2 + P3, Si = {pj + Pk)'^ and qi = pj — Pk, where {i,j, k) represent 
cyclic permutations of (1,2,3). The factor T'^'^ denotes the expression g'^'^ — a^oy ja}. The factors 
Bp{si) denote Breit-Wigners describing the corresponding p -^ tttt^ amplitudes. Finally, we have 
included an additional form factor Fr^^, which represents the effect of the finite size of the ai 
meson in its decay to pn. We take this form factor to have the form 



FR^ih 






(5) 



where hi is the momentum of the decay products, the p and the tt in this case, in the a\ rest frame. 
The parameter Ri is proportional (by a factor of \/67ic, see Ref. |^) to the root mean square (r.m.s.) 
radius of the a\. We note that expressions for jf must be symmetric with respect to interchange 
of 'n:\ and t^\ since these are indistinguishable. 

In our analysis of substructure in the r — > Vj-'in decay, we consider the following amplitudes: 



f. 


s-wave amplitude for 1+ 


f2 


s-wave amplitude for 1+ 


f. 


d-wave amplitude for 1+ 


3l 


d-wave amplitude for 1+ 


i^ 

h 


p-wave amplitude for 1+ 


3', 


p-wave amplitude for 1+ 


fl 


p-wave amplitude for l"*" 



p(770)7r (denoted as pvr), 

p(1450)7r (denoted as /o'vr), 

/o(770)7r, 

/9(1450)7r, 

/2(1270)7r, 

/o(400 - 1200)7r (denoted as o-vr) 

/o(1370)7r. 



(6) 



An explicit parametrization of the amplitudes j^ is given by Eqn. ^^. With these definitions, the 
constants /3j have dimensions of [GeV]^, where the exponent depends on the amplitude. In our fits 
we specify /3i = 1, such that the couplings for the other amplitudes are determined relative to the 
s-wave yOTT coupling. The parameters used to describe the resonances appearing in the jf above are 
given in Table |, while the Breit-Wigner form used here is given in Appendix ^ by Eqn. A7. 
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TABLE I. 


Resonance parameters for the intermediate 


states 


as used 


in 


the substructure fits. 


Y 


[GeV] 




[GeV] 




Ref. 



p(770) 0.774 0.149 p| 



p(1450) 1.370 0.386 p| 



/2(1270) 1.275 0.185 O 



a 0.860 0.880 |0| 



/o(1370) 1.186 0.350 [20| 



In Eqns. Q and A3, we have constructed Lorentz- invariant amphtudes so as to make contact 
with the resonant components of the substructure. In contrast with a formulation based on angu- 
lar momentum eigen functions, these amplitudes are only approximately associated with a specific 
angular momentum quantum number L, and hence we have employed lower case letters to iden- 
tify the primary value of L. Thus, for example the lowest dimensional Born amplitude for /97r, 
the Lorentz-invariant s-wave amplitude, contains a small D-wave component (see for example, 
Refs. i§). 

The selection of the amplitudes ji , . . . , jj is in part based on experience gained in early attempts 
to fit the data. It is also in part motivated by the unitarized quark model of Tornqvist [EG]. The 



resonance parameters of the broader a and /o(1370) mesons are taken from application of this 



model to existing data [20|. We have also performed fits with additional amplitudes, namely the 



axial vector /o(980)7r and pseudoscalar ir' -^ pn and an. These are discussed in Section IVD 



B. Model for the Stt mass spectrum 

The conventional understanding of r — > j^,- Sit decay is that it proceeds through creation of the 
lowest lying axial vector meson, the ai(1260). Since radial excitations may also be present, we 
replace the Breit-Wigner function Ba^ (s) appearing in Eqn. ^ by a modified function that includes 
a possible a[ admixture. 

B{s) = Ba,{s) + K ■ B,,{s) = YTT—- ^^^^^rT-\ 

+ 2 — ^<7T' (') 

where k is an unknown complex coefficient. The a'^ meson is predicted in a flux-tube-breaking 



model [21,22[ with a mass of m-oa' = 1-820 GeV. Experimental indications [23[ suggest a mass of 
1.7 GeV and a width of 0.3 GeV. One impact of introducing the a'l in this way is that the coupling 
constants /3i in Eqn. ^ will necessarily vary with s. We will return to this issue later in this article. 
In Eqn. 1^, the function w.^^(s) is the running mass P,pO[, 



m 



iM=<a,+^\s), (8) 



where (5^(s) is the mass shift function, 
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The mass shift function is renormahzed such that 

"lai(s)L=m2 =moai- (10) 

ai 

The bare mass rrioai is chosen to be the resonance mass by requiring that the total width r"^j(s) 
at s = rriQ ^^ is equal to the nominal width Fq ai • 

The \/i-dependent behavior of the ai width, and consequently of its mass, requires knowledge 
of the underlying substructure, not just for ai — > Sir, but also for decays to other channels such 
as ai — > KKtt [via ai -^ K*K and oi -^ /o(980)7r]. Considering only these contributions, the oi 
width can be written as 

where 50^(3,^) denotes the coupling of the oi meson to the Svr system, the Jai{x) denote the relative 
coupling of X to the oi meson, and the F^^ denote the reduced widths. 

The ai — > Svr partial width can be expressed in terms of the amplitudes for the hadronic current 
jf , as defined by Eqn. |3[ Specifically, the reduced widths Tg^o^T ^^^ '^2n^Tr± ^^^ 






where d^sj^ denotes 3-body phase space for the a^ — > vr^vr'^vr'^ decay. We determine f 2^0,^^ ^^' 
merically using the output from the substructure analysis in which the jf are specified according 
to Eqn. |A3| in Appendix lAl Similarly, for r2^^^± we make use of isospin symmetry to infer the 
corresponding all-charged amplitudes from our analysis of the 7r~7r''7r'^ substructure. 

The K*K and /o(980)7r partial widths contribute as thresholds in the y^-dependence of the 
ai width. We determined these from expressions for the s- and p-wave amplitudes, respectively, 
making use of the narrow width approximation for the /o and K* . The relative couplings for these 
contributions are left as free parameters to be determined from the data, along with the ai pole 
mass and width. 

IV. ANALYSIS OF DALITZ PLOT AND ANGULAR VARIABLES 

The primary goal of this analysis is to characterize the contributions to the substructure of the 
oi — > Svr decay, as well as the parameters describing ai line shape itself, including the question 
of possible radial excitations. Two separate analyses are carried out to address these two issues, 
however it is important to realize that they are closely coupled. 

First, the integration over the Dalitz plot needed to specify the mass-dependence of the ai width 
as well as the running of the oi mass requires the amplitudes participating in the Svr hadronic 
current to have been determined. On the other hand, the question of whether an a'l resonance 
also contributes to the Svr mass spectrum affects the way one would choose to parametrize the 
substructure. 

Practically, it is not feasible to fit the Svr mass spectrum and the hadronic substructure simulta- 
neously. We have elected to determine first the substructure in a way that is mostly independent of 
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the Svr mass spectrum. Then in a second step, using the results on f 3^(s) obtained in the substruc- 
ture fits, we measure the ai resonance parameters from the Stt mass spectrum. In this section we 



describe the substructure fits, while the fits to the Svr mass spectrum are described in Section VI. 



A. Fitting method 



To determine the contributions to the substructure in r ^ ^'t-Stt, we perform unbinned maxi- 
mum likelihood fits using as input the measured 3-momenta of the three pions in the decay, and 
the energy of the decaying r lepton. The latter is known to be the beam energy in the absence of 
radiative effects. With knowledge of the particle masses (we take the mass of the r neutrino to be 
zero), these inputs comprise a complete kinematical description of the decay, with the exception 
of: (1) deviations in the r energy due to initial state radiation (ISR), (2) the azimuthal orientation 
if^'^ of the r flight direction relative to the measured momentum vector of the Svr system, and (3) 
smearing due to scattering and detector resolution, the beam energy spread, and radiative effects 
other than ISR. 

Following the discussions in Sec. Ill and Appendix ^, and ignoring the sources of smearing 
described in item (3) above, we construct the likelihood function. The numerator da{pi,p2,P3) of 
the likelihood £ = da/ J da is 



da 



P (cos 



?3;',COSI 



e'^^Eif, \p^\), ipl^) X f{p^) X (S^, + iK^A^,)}d%d^l^ 



Y.if^di'm jT) X \Bis) 



dcosO'^l'dift'dEifd^s^ds, 



(13) 



where we integrate over the unmeasured Ur [(p^ above) and ISR photon {jp^) degrees of freedom, 
such that da is a function of the measured degrees of freedom. For illustrative purposes we represent 
these in the above by the squared 37r invariant mass (s), the energy and orientation of the 37r system 
in the laboratory {E^^, 9^^^, ip'f^), and the 3-body phase space element (^$37^). The phase space 
factor can be expressed in terms of the Dalitz plot variables si = M^_ and S2 = M^_ and the 
Euler angles representing the orientation of the Sir decay plane in the 37r rest frame. The symbols 
S^i, and Afj^iy represent the symmetric and antisymmetric terms in the lepton tensor. The factor 
fiP'j) denotes the factorized ISR photon probability distribution. Finally, we have also included the 
r-pair production dynamics, the effect of which is to make non-uniform the probability distribution, 
denoted by the factor P, for the azimuthal angle (p^'^ . The polar angle 6^"^ between the r direction 
and the Sir system appearing in this factor is determined by -Ebearm E^l^ and \p^\. The r neutrino 
helicity hiy^ and the complex coupling constants Pi of the hadronic amplitudes are the fit parameters. 
The above integral is computed using a reverse Monte Carlo technique |^,24|. In this method, 
for each event in the data we generate a sample of trial MC events which are designed to have 
precisely the measured values for the pion momenta, but which have unmeasured quantities de- 
termined randomly according to the factorized distributions for ISR photons and the unknown 
azimuthal angle (f^^ . The integration is performed using trial events that possess internally consis- 
tent kinematics. We remove data events for which the number Uhu of these successful trials is low 
so as to maintain high precision on the integration. This requirement also tends to preferentially 
remove background events. 
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TABLE II. Background contributions (%) to the lepton-tagged sample in slices of the 3tt in- 
variant mass ^/s after the cuts on y^ > 0.6 GeV and n^it- 



fake it' 







v-rT^ 3tt 







Ut-K vr vr 



0^0 



T — > Vt-KsT^ 



Bin 1 


0.6-0.9 


GeV 


17.2 ±1.5 


14.0 ±1.6 


2.1 ±0.6 


8.7 ±1.3 


Bin 2 


0.9-1.0 


GeV 


12.9 ±0.8 


5.5 ±0.6 


0.7 ±0.2 


1.5 ±0.3 


Bin 3 


1.0-1.1 


GeV 


9.5 ±0.5 


4.3 ±0.4 


0.5 ±0.1 


0.1 ±0.1 


Bin 4 


1.1-1.2 


GeV 


6.8 ±0.4 


3.0 ±0.3 


0.6 ±0.1 


0.1 ±0.1 


Bin 5 


1.2-1.3 


GeV 


7.1 ±0.5 


2.3 ±0.3 


0.7 ±0.2 


0.0 


Bin 6 


1.3-1.4 


GeV 


6.5 ±0.6 


1.8 ±0.3 


0.4 ±0.2 


0.0 


Bin 7 


1.4-1.5 


GeV 


6.6 ±0.9 


0.9 ±0.3 


0.2 ±0.2 


0.0 


Bin 8 


1.5-1.8 


GeV 


6.8 ±1.3 


0.3 ±0.3 


0.2 ±0.2 


0.0 


Bin 1-8 


0.6-1.8 


GeV 


8.5 ±0.2 


3.5 ±0.2 


0.6 ±0.1 


0.5 ±0.1 



To be insensitive to details of the 37r mass spectrum (addressed in Section [V^ ), we subdivide 
the data in fine bins (25 MeV) of ^/s and calculate the normalization N = f da of the likelihood 
separately for each bin j: 



Ni 



'j+i 



e • da 
dcos ei^f d^f^ dEl^f d$3^ ds 



dcOSe',fd^l^^'d4fd<l>3n 



ds , 



(14) 



where e denotes the detector efficiency. Over the bin width As = Sj+i — Sj, \B{s)\ is approximated 
to be constant, and thus cancels in the likelihood. The normalization integrals are computed using 
factorization-based Monte Carlo events that have been passed through the full detector simulation 
as described in Section II C| . 



B. Treatment of backgrounds 



In addition to the likelihood for signal events defined by Eqn. |T^, we also include the four main 
background sources listed in Table ||. There, the background fractions, estimated from the r MC 
sample for the i/47r, uKittt and vKstt modes, are tabulated in slices of ^/s so as to illustrate the 
dependence. 

-> Vt-p~ events where a spurious vr'^ has been recon- 



Events with fake tt'^'s tend to be r 



tructed from clusters associated with radiative photons, shower fragments from the interaction of 
the charged vr in the detector, or other accidental activity in the calorimeter. The likelihood dis- 
tribution for the fake tt^ background is approximated from data by the Dalitz plot distribution of 
events populating the vr" mass side bands. 

For the r~ -^ ^'^7r~37r'' background the reverse Monte Carlo procedure is modified to simulate 
a lost vr''. The 47r matrix element is not well measured. We consider models in which the 47r system 
arises via the p(1450) resonance, where we simulate either (1) p'~ -^ a^-K^ (s-wave) -^ p^ir^-K^, or 
(2) p'~ — > p~a (s-wave), or a combination thereof. The Dalitz plot projections from these models 
are very similar. In addition, the goodness of fit varies little with the choice of model. In the fits 
reported here, we used the p'~ — > p~a model. 
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The background t~ -^ Vt-K~'k^t:^ is modeled by the decay chain r~ — > u^-Ki , K^ -^ 
K^^ir^ (s-wave), where the Ki meson is parametrized by a superposition of the i^i(1270) and 
i^i(1400) Breit-Wigner functions. Finally, the r~ -^ j/^i^^vr" background is parametrized by the 
decay chain t~ — > v-j-K*~ , K*~ — > K^t^~ (p-wave). The mass distribution for the Ks -^ tt'^tt'^ 
decay is parametrized by a Gaussian, where the mean and the width are taken from data. 

With the inclusion of these backgrounds the likelihood function is: 

C = {I — af^g — a^T, — aKiTTT — a Ks-k)^ signal 

+ Oif^o^f^o + 0!.4:n^A7T + OiKnir^Kmr + OiKsn^Ks'TT (15) 

The background fractions Oj depend on ^/s and are taken from Table |l|. 

C. Results 

In this section, we report on the fits to the substructure in t~ — > UrT^^-K^TT^ decays. Given the 
complexity of the fitting procedure, we use only the lepton-tagged sample since the backgrounds 
from multihadronic events and other r decays are smaller, particularly in the high Svr mass region. 
We have performed many fits, including various amplitudes and employing differing assumptions. 



Here, we present results from one fit based on the model described in Section III A, with certain 



parameters fixed as described below. Results obtained when these parameters were varied are given 



in Section IV D 



The resonances Y shown in Table | are implemented in the fit in amplitudes for A — > Ytt, 
where A represents an axial vector system. As mentioned above, we compute the normalization 
of the likelihood function in bins of y/s so as to be insensitive to the resonant content of the Stt 
system. In addition, the couplings (3i could vary as a function of ^/s. This could be the case if, for 
example, several resonances contribute to the Svr system. In our nominal fit, we constrain the /?, to 
be independent of ^/s. For simplicity, we also consider the l"*" system to be point-like, i.e., we set 
-Rj = in Eqn. |5|, with the result that Fr. = 1. Finally, we fix the Ut helicity h,^^ to its Standard 
Model value of —1. Thus, our fit parameters consist of twelve real numbers: the moduli \(3i\ and 
phases (pp. of the couplings, for i = 2-7. Fits with h^^ floating are discussed in Section ^. 



The results from this nominal fit are summarized in Table |I|. The measured likelihood is 
224259, while that expected is 224340 it 214. The difference, — 0.4 cr, indicates an acceptable good- 
ness of fit. As a function of ^/s, it is (in units of standard deviations): +2.2, —1.3, —0.3, —1.5, +0.9, 
—0.7, —0.6 and +1.2 in the eight slices of ^/s defined in Table ||. The significance of each amplitude 
is determined by repeating the fit with that amplitude excluded. Dalitz plot projections from the 
fit are shown in Figs. |^ and ^ in slices of ^/s, overlaid with the corresponding data distributions. 
The Dalitz plots themselves are shown in Fig. g. A discussion of the results follows in Section IV F . 



For now, we note the large contributions from channels involving isoscalars, in particular air with 
a significance of 8.2a. 
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TABLE III. Results of the nominal fit for the moduli |/3j| and phases (/>^. of the coefficients for 
the amplitudes listed in Eqn. |6|. The two errors shown are statistical and systematic respectively. 
The branching fractions B are derived from the squared amplitudes (using the values of |/3i|), and 
are normalized to the total t~ -^ UrTr'ir^ir^ rate. These do not sum to 100%, due to interference 
between the amplitudes. 



Signif. 



m 



(t>pjT^ 



B fraction(%) 



p 


s-wave 


— 


p(1450) 


s-wave 


1.4a 


P 


d-wave 


5.0a 


/9(1450) 


d-wave 


3.1a 


/2(1270) 


p-wave 


4.2a 


a 


p-wave 


8.2a 


/o(1370) 


p-wave 


5.4a 



1 

0.12 ±0.09 ±0.03 
0.37 ±0.09 ±0.03 
0.87 ±0.29 ±0.06 
0.71 ±0.16 ±0.05 
2.10 ±0.27 ±0.09 
0.77 ±0.14 ±0.05 




0.99 ±0.25 ±0.04 

-0.15 ±0.10 ±0.03 
0.53 ±0.16 ±0.06 
0.56 ±0.10 ±0.03 
0.23 ±0.03 ±0.02 

-0.54 ±0.06 ±0.02 



68.11 
0.30 ±0.64 ±0.17 
0.36 ±0.17 ±0.06 
0.43 ±0.28 ±0.06 
0.14 ±0.06 ±0.02 
16.18 ±3.85 ±1.28 
4.29 ±2.29 ±0.73 
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3100199-003 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

2 
s or s (GeV ) 

FIG. 3. Dalitz plot projections: distributions in squared tt^tt^ mass si and S2 (two entries per 
event). The data are represented by the filled points. The solid line is the fit result. Background is 
represented by the shaded histograms. The lightest shaded histogram is the sum of the backgrounds, 
while the darker histograms show the backgrounds separately. Plots (a) to (h) correspond to slices 
in ^= 0.6-0.9, 0.9-1.0, 1.0-1.1, 1.1-1.2, 1.2-1.3, 1.3-1.4, 1.4-1.5, 1.5-1.8 GeV. 
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3100199-004 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

Sg (GeV^) 

FIG. 4. Dalitz plot projections: distributions in squared vr'^vr'^ mass S3 (one entry per event). 
The data are represented by the fihed points. The soUd hne is the fit result. Background is 
represented by the shaded histograms. The lightest shaded histogram is the sum of the backgrounds, 
while the darker histograms show the backgrounds separately. Plots (a) to (h) correspond to slices 
in ^= 0.6-0.9, 0.9-1.0, 1.0-1.1, 1.1-1.2, 1.2-1.3, 1.3-1.4, 1.4-1.5, 1.5-1.8 GeV. 
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Bin 2, Fit Result 



0.2 0.4 0.6 0.2 0.4 0.6 0.2 0.4 0.6 0.2 0.4 0.6 



CM 



1.0 
0.8 
0.6 
0.4 
r- 0.2 



0) 

O 

(A 1.0 



T—' — I — ' — \ — ' — I — ' 

Bin 3, Data 



— 1 1 1 1 1 1 1 1 r— 

Bin 3, Fit Result 



J I I I L 



0.2 0.6 1.0 0.2 0.6 1.0 0.2 0.6 1.0 0.2 0.6 1.0 



-' — I — ' — I — ' — \ — ' — r 

Bin 4, Data 



J I I I L 



— 1 1 1 1 1 1 1 1 r— 

Bin 4, Fit Result 



.£: 



J I I I I , L 



0.5 



Bin 5, Data 



Bin 5, Fit Result 



Bin 6, Data 



_!_ 



Bin 6, Fit Result 



_!_ 



0.5 1.0 



0.5 1.0 



1.5 

1.0 

0.5 





Bin 7, Data 



Bin 7, Fit Result 



0.5 1.0 

I ' ' ' ' I ' ' ' ' I 
Bin 8, Data 



0.5 1.0 1.5 

I ' ' ' ' I ' ' ' ' I ' ' ' 
Bins, Fit Result 



0.5 1.0 1.5 0.5 1.0 1.5 0.5 1.0 1.5 0.5 1.0 1.5 2.0 

s^ (GeV^) 

FIG. 5. Dalitz plot distributions for data and fit result. Here si is taken to be the higher 
of the two values of M^_ o in each event. Bins 1 through 8 correspond to slices in a/s = 0.6-0.9, 
0.9-1.0, 1.0-1.1, 1.1-1.2, 1.2-1.3, 1.3-1.4, 1.4-1.5, 1.5-1.8 GeV. 
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D. Modifications to the default model 



With the model described in Section III A, we have obtained a good fit to the Dalitz plot 



distributions. In this section we describe fits to variations of the model. 



1. Uniformity of amplitude coefficients across ^/s 

The assumption that the coefficients Pi for the various substructure amplitudes are independent 
of ^/s may not be correct. They would not be constant if, for example, more than one Stt resonance 
were present. We have performed fits allowing the f3i to vary; the results from one such fit are 
plotted in Fig. q. In this fit, we considered fewer amplitudes so as to limit the number of fit 
parameters. They are pir (d-wave), /2'7r (p-wave) and air (p-wave), in addition to the dominant s- 
wave pTT contribution. Also, for this fit, we take the p resonance to be the sum of p{770) and yo(1450) 
amplitudes, with the p(1450) admixture fixed according to studies of r~ —>■ i/t-tt^tt'' decay |l^. The 
goodness of fit is acceptable: the measured likelihood minus that expected is +0.2a. 

The behavior of the moduli \Pi\ are consistent with uniformity across ^/s. However, we note that 
the signficance of the /27r contribution is greatest in the highest mass slice. We also see elevated 
contributions from the cnr and d-wave pn amplitudes in the high- mass slices, although these are 
not statistically significant. 



3100199-006 




0.75 1.25 1.75 0.75 1.25 1.75 

S''^ GeV 

FIG. 6. Results for the moduli and phases of coefficients Pi of amplitudes from a fit in which 
these are allowed to vary as a function of the Svr invariant mass. In this fit, only three amplitudes 
are considered in addition to the s-wave pir amplitude: pir (d-wave), /27r (p-wave) and air (p-wave). 
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2. Importance of isoscalar contributions 



This analysis is the first study of the ai in r decay to consider contributions from scalar 7 = 
mesons [a and /o(1370)]. In addition, our fits return a significant /2(1270)7r component. Although 
these channels are expected to be present, a demonstration of the validity of the fit results is 
desirable given the complexity of both the model and the fit procedure. To help visualize their 
collective importance in describing the substructure, we have performed fits excluding the three 
amplitudes js, jg, and jV that involve isoscalars. The Dalitz plot projections from one fit performed 
in this way are presented in Appendix y for comparison with those from the nominal fit. We 
comment further on the impact of the large isoscalar contributions in Section IV F|. 



3. Consideration of finite size of the ai meson 



In our nominal fit, we set the ai radius to zero, such that the form factors Fr. in Eqn. g are 
just unity. We find that good fits can also be obtained with non-zero values for Ri. In Fig. ^, we 
plot the differences in the minus log likelihood values from fits in which all Ri are set to some value 
R. The best fit is obtained with R = 1.4GeV~ . We present the results from this fit in Table VII 
in Appendix |C|. We will return to this issue in the context of the fits to the Svr mass spectrum in 
Section 0. 




R(GeV '*) 



FIG. 7. Dependence of —2 In C on the value of the meson radius parameter Ri 
of the nominal {Ri = 0) substructure fit. 



R in variations 



4- Inclusion of ^ Svr amplitudes 

We have also performed fits including pseudoscalar contributions, namely vr' -^ pn (p-wave) 
and vr' -^ an (s-wave). These amplitudes will necessarily have a different i/s-dependence from 
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those associated with axial vector production. To account for this, we assume a Breit-Wigner form 
for the tt' with a constant mass and width of 1.300 and 0.400 GeV respectively, and use the results 



from one of the fits presented in Section VI for the ai parameters. 

In fits with each of these amplitudes included separately we find no statistically significant 
contributions. We obtain the following 90% CL limits: 

B{t -^ vtt' -^ upiT -^ u3tt) < 1.0 X 10"^, (16) 

B{t -^ vTi' -^ van -^ z^37r) < 1.9 x 10"^. (17) 



5. Inclusion of other 1+ -^ Svr amplitudes 

In addition to the axial vector amplitudes if , • • • , jy , we performed fits including a contribution 
from /o(980)7r (p-wave). None of these fits returned a signficant coupling for this amplitude. Further 



discussion of possible /o(980) contributions appears in Section VIC 4 



6. Variation of a meson resonance parameters 

By virtue of its low mass and large width, there is much uncertainty regarding the resonant 
shape of the a meson. For simplicity, we have elected to characterize it using a Breit-Wigner form. 



with its mass and width taken from the model of Tornqvist [20|. We have not considered alternative 
forms, or explored extensively the range of possible resonance parameters. However, in view of the 
large a contribution observed in this analysis, we have attempted to ascertain whether our data 
are sensitive to variation of its properties. 

We have refitted the data with a range of input values for the mass and width of the a. Of the 
values we considered, the best fit was obtained with moo- = 555 MeV and Foo- = 540 MeV. The 
value of — 21n£ for this fit was 224216. This is 43 units below that for the nominal fit, but is still 
consistent with expectations given the statistics of the data sample. Using the smaller values of 
rriQfj and Fqo- has an impact on the values of f3i obtained. The main trend is a relative change of 
20-40% in the ai branching fractions, which are smaller for the cttt and /o(1370)7r channels and are 
larger for the p'vr and /2(1270)7r channels. 

E. Systematic errors 



The systematic errors shown in Table III are based on estimates of the uncertainties arising from 
the following sources: Monte Carlo statistics, background fractions and modeling, dependence of the 
acceptance on the kinematical observables used in the fit, and detector resolution. The uncertainties 



due to these sources are given in Table IV 



The error due to Monte Carlo statistics is based on the variance of results obtained from six 
separate fits, each using one sixth of the Monte Carlo sample for the normalization of the likelihood 
function. Fits performed after varying the background fractions and model (in the case of the 
Vt 7r~ Stt*^ channel) within reasonable limits were used to estimate the error associated with this 
source. To estimate the error associated with acceptance, the Monte Carlo was used to parameterize 
acceptance as a function of the charged and neutral pion momenta as well as the opening angles 
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TABLE IV. Systematic errors on hadronic substructure for the nominal fit. 







/)(1450) 


P 


p(1450) 


/2(1270) 


a 


/o(1370) 






s-wave 


d-wave 


d-wave 


p-wave 


p-wave 


p-wave 


Monte Carlo 
statistics 


A{B fraction) (%) 


±0.162 


±0.058 


±0.045 


±0.017 


±1.201 


±0.701 


A(IAI) 

A(0ft/vr) 


±0.029 
±0.038 


±0.028 
±0.027 


±0.044 
±0.058 


±0.047 
±0.026 


±0.081 
±0.011 


±0.049 
±0.020 




A{B fraction) (%) 


±0.027 


±0.008 


±0.042 


±0.014 


±0.391 


±0.157 


background 


A(IAI) 


±0.006 


±0.005 


±0.037 


±0.010 


±0.036 


±0.018 




A(<Aft/vr) 


±0.005 


±0.007 


±0.018 


±0.009 


±0.007 


±0.006 




A{13 fraction) (%) 


±0.024 


±0.014 


±0.006 


±0.002 


±0.178 


±0.097 


efficiency 


A(IAI) 


±0.004 


±0.007 


±0.008 


±0.002 


±0.016 


±0.007 




A(<Aft/vr) 


±0.013 


±0.002 


±0.009 


±0.001 


±0.006 


±0.008 


detector 
resolution 


A(^ fraction) (%) 


±0.018 


±0.006 


±0.004 


±0.002 


±0.083 


±0.065 


A(IAI) 


±0.001 


±0.002 


±0.007 


±0.001 


±0.011 


±0.006 


A(</>ft/^) 


±0.002 


±0.004 


±0.006 


±0.002 


±0.004 


±0.005 




A (5 fraction) (%) 


±0.17 


±0.06 


±0.06 


±0.02 


±1.28 


±0.73 


total 


A(IAI) 


±0.03 


±0.03 


±0.06 


±0.05 


±0.09 


±0.05 




A(<Aft/vr) 


±0.04 


±0.03 


±0.06 


±0.03 


±0.02 


±0.02 



between these particles. Reasonable deviations from these parameterizations were used to reweight 
events entering the fit, and the resulting variations in fit parameters were taken as the systematic 



errors. Finally, the likelihood function in Eqn. 13 does not take into account resolution effects. The 
effects of modifying it to include resolution smearing based on errors in track parameters for the it~ 
and in photon energies and directions for the n^'s was used to estimate the error from this source. 
For all fit parameters, the error due to limited Monte Carlo statistics dominates the systematic 
error. 



The results given in Table III are meaningful only in the context of the model used to 



parametrize the substructure. Different models yield results that can differ significantly from our 
nominal fit results. Given this plus the unfeasibility of examining all possible models, we have not 
attempted to assign a systematic error associated with model dependence. 

F. Discussion 

The results of the fits for the substructure can be summarized as follows: 

• The pvr s-wave amplitude with a branching fraction of around 70% is dominant, as expected. 

• With the exception of the yo(1450)7r s-wave amplitude, all amplitudes included in the nominal 
fit contribute significantly to the Sir hadronic current. In other fits, we find no evidence for 
contributions from ai — > /o(980)7r, or from t~ — > z^t-vt'". 

• The isoscalar mesons /2, /o(1370), and a contribute with a combined branching fraction of 
approximately 20% to the 37r hadronic width. In particular, the a meson with a significance 
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• 



of ^ 8a cannot be neglected. It shows up strongly as part of the broad enhancement at the 
low end in the S3 (M^o 0) distribution in the slices of ^/s shown in Fig. ^(g) and (h), but is 
significant in all slices. 

The p{14:50)tt state shows up more strongly in the d-wave amplitude than in the s-wave 
amplitude. 



The last point above may have implications regarding a possible a'l contribution. One expects 
that an a'l contribution induces -v/s-dependent couplings. In the fits allowing Pi to vary with ^/s, 
we found that (1) the goodness of fit is not significantly better, and (2) the values of /3j(s) are 
roughly consistent with being constant. On the other hand, according to the flux-tube-breaking 
model of Refs. [2l|,22|, the a'^ meson prefers to decay to pir by D-wave rather than S'-wave, and 



the /9(1450) is preferred over the p{770). Thus, it is possible that the measured p{1450)tt d-wave 
amplitude could be induced by an a'^. The suggestions of enhanced /2vr and an contributions at 
large y/s are also consistent with the hypothesis of an a'l. However, the statistics of the present 
data sample are not sufficient to resolve this question with the substructure fits. 

As a test of our fit results, we have compared the Dalitz plot distributions from a sample of 
background-subtracted t~ — > ^'T-vr^Tr^vr^ events with the isospin prediction based on the results 
from the nominal fit to the Tr^vr^vr'^ mode. The backgrounds were estimated from the generic r 
Monte Carlo sample. The dominant backgound t~ — > ^'^7r~7r~7r+7r'' is simulated in this sample 
with the model implemented in TAUOLA, containing Tr~uJ as well as [pTTir]^ (in various charge com- 
binations) substructure. The Dalitz plot projections are shown in Fig. p. The observation that the 
hadronic current for the all-charged mode is well described by our results for the Tr^vr^vr'^ mode 
provides a critical corroboration of our measurements. This is particularly important for the am- 
plitudes involving isoscalars since they enter the all-charged mode with the opposite sign relative 
to the other amplitudes. A full analysis of the high-statistics all-charged mode is under way and 
will be presented in the future. 

The world average values IQ for the r~ — > z^T-vr^Tr+vr^ and r" -^ UrTr'n^n^ branching frac- 
tions are (9.23 it 0.11) % and (9.15 it 0.15) %, respectively. Their near equality is consistent with 
expectations from isospin symmetry, if the decays were to proceed exclusively via pir or p'n. One 
consequence of the presence of isoscalars in oi decay is the possibility of upsetting this expectation. 
However, due to interference the two modes contribute nearly equally to Tl^f{s) (see Fig. ^(b) in Sec- 



tion VIB). The divergence of these contributions at high values of y/s is damped in the decay rate by 
the falling of the ai line shape, as well as by phase space and weak interaction dynamics in the r de- 
cay. Furthermore, the residual preference for the all-charged mode at high ^/s is compensated for by 
the larger phase space available for the 7r~27r^ mode at low y^. Quantitatively, the ratio of branch- 
ing fractions is predicted from this analysis to be B{t^ -^ Vt^^ tt^ n^) / B{t^ -^ Ut-tt^it^tt^) = 0.985, 
in agreement with the ratio 1.009 ib 0.020 obtained from the direct measurements. 



Finally, the branching fractions reported in Table [H are the t~ decay branching fractions 
relative to the total i't'k~'k^'k^ rate. These differ from the ai branching fractions due to the 
weighting of the oi line shape by factors associated with r^ — > Vra^ weak decay. The oi branching 
ratios ;S(af -^ \^''^]^ -^ tt^ti^tt^) /B{ai —>■ tt'tt'^tt^) are given in Table ^ 
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s^ or s^ (GeV ) 

FIG. 8. Background-subtracted squared vr^vr^ mass (si and S2) spectrum for the three charged 
pion mode (two entries per event). The data are represented by the fihed points. The sohd hue is 
the isospin prediction based on our fit to the 7r~7r'^7r'^ mode. 



V. DETERMINATION OF THE SIGNED u^ HELICITY 



In the fits reported in the previous section, the r neutrino (anti-neutrino) heUcity was fixed to 
the Standard Model value of h^^ = — 1 {h-p^ = -|-1). However as first pointed out by Kiihn and 
Wagner [^, interference between the two 7r~TT^ systems gives rise to a parity violating term in the 
squared matrix element for the decay r~ -^ Vt-ti~'k^tt^. This permits determination of the sign, 
as well as the magnitude, of the neutrino helicity. Including h^^ as an additional free parameter 
to the nominal fit described in the previous section, and assuming invariance under the combined 
charge-conjugation and parity {CP) operation, we find h^^ = —1.00 it 0.12. In this fit, the values 
of (3i are affected at a negligible level. Including only the s-wave pvr amplitude in the model for the 
substructure yields a poor fit, with h,y^ = —0.73 ± 0.09. 

To investigate the model dependence entering this measurement we also performed fits for 
substructure and hy^ with different input parameters. While our nominal fit was performed with 
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TABLE V. Branching ratios for a^ decay into vr vr'^vr'^ via intermediate states shown, relative 

to the total af -^ tt^tt^tt^ rate. The errors shown are statistical only. 

Amplitude Branching ratio (%) 

pTT s-wave 60.19 

/9(1450)7r s-wave 0.56 ± 0.84 

pn d-wave 1.30 ib 0.60 

p{U50)tt d-wave 2.04 ± 1.20 

/2(1270)7r p-wave 1.19 ±0.49 

crvr p-wave 18.76 ± 4.29 

/o(1370)7r p-wave 7.40 ± 2.71 



the oi radius set to zero, non-zero values also gave good fits, as noted in the previous section. Using 
the best fit value of Ri = 1.4 GeV~ , we obtain h^^ = —1.03 it 0.13 . As a best estimate of h^^ given 
the dependence on input assumptions, we average this value with the Ri = result to obtain 

h^^ = -1.02 ± 0.13 (stat.) ± 0.01 (syst.) ± 0.03 (model) , (18) 

where the uncertainty due to the model dependence is estimated by the difference between the values 
from these two fits. This result agrees with other determinations |jl| of the sign and magnitude of 
hu^, as well as with the Standard Model value of —1. 

The systematic error given for h^^ was determined in a fashion similar to those in the sub- 
structure analysis. The sources contributing to this error are: Monte Carlo statistics (±0.005), 
background determination (±0.010), dependence of the acceptance on the kinematic observables 
used in the fit (±0.003), and detector resolution (±0.004). 

We have also looked for possible CP non-conservation by determining h^,^ and h^^ separately. 
Defining a CP-violating asymmetry 

^CP = i^itl^^l , (19) 

we find Acp = —0.08 ± 0.13, where the error is dominantly due to statistics. 



VI. ANALYSIS OF THE 7r-7r°7r° MASS SPECTRUM 

In this section, we describe fits to the 3tt mass spectrum performed to extract resonance param- 
eters of the ai meson. The results from the substructure fits are used as inputs for determination 
of the running of the oi mass and width. 



A. Fitting method and assumptions 

The ai resonance parameters are determined from a x^ fit to the background-subtracted and 
efficiency-corrected vr^vr'^vr'^ mass spectrum. The fit is performed over a range from M^t^ = 0.600 
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to 1.725 GeV in bins of width 25 MeV. As in the substructure analysis, we perform a variety of 
fits, reflecting different models and assumptions. We perform fits to the all-tagged sample, as well 
as to the lepton-tagged subsample, to benefit from the higher statistics. 



For our nominal fit, we specify the following version of the model described in Section |IIIB 
(1) turn off a'l contribution, i.e., set k = in Eqn. ^; (2) turn off form factors describing finite 
size of the oi, i.e., set Ri = 0; (3) include the K*K threshold, but not the /o(980)7r threshold, in 
determining Tll^[s)] and (4) assume the running ai mass to be flat as a function of ^/s. We have 
also performed fits in which various of these specifications are modified, and obtain satisfactory 
results under a variety of configurations. With the above specifications, the fit contains three free 
parameters in addition to an overall normalization: the pole ai mass moai , the ai — > Svr coupling 
9ai(37r)) and the relative ai — > K*K coupling ^ai{K*K)- We derive from these parameters the pole 



oi width Foai using Eqn. 11 



Some comments on the above choices are in order. The -v/i-dependence of the oi total width 
depends strongly on assumptions. Inclusion of the K*K channel is motivated by observation of the 
decay t~ -^ Vr[K*K]^ [Q], although it is not well-determined as to how much of this comes through 
the axial vector (rather than vector) weak current. As we have no evidence for the /o(980)7r channel 
in the substructure fits, we have omitted its possible contribution in our nominal fit here. However, 
this and other thresholds may be present. For example, a possible ai — > /i(1285)7r channel, as 
suggested by the recent obervation of this system in r decay [^, would open up near the K*K 
threshold. Thus, the value for F^*^ returned from our fit can not be strictly interpreted as just 
the ai -^ K*K partial width. 

The running of the oi mass is even more problematic since the upper limit of integration (over 
ds') in Eqn. ^ is infinity, and thus the integral will include effects from channels that open above 
the T mass, and are therefore not directly measureable in r decay. Furthermore, a damping of 
the amplitudes, such as that provided by the form factors Fr. , is needed so that the integral can 
converge. As a result, allowing the ai mass to run is practical only in models where the Ri are non- 
zero. In such models, the effect of additional thresholds at high ^/s is to flatten the y^-dependence 
of the running mass. Thus, we expect the running mass to be closer to a constant than we would 
predict from Eqn. |8| with known thresholds. Setting Ri = and taking a constant ai mass may 
not be rigorous, however the resulting model is simplified. 

B. Results 

The results obtained from our nominal fit to the Svr resonance shape parameters are shown in 
the second column of Table ^. The x^ for this fit is 39.3 for 41 degrees of freedom. Fits to just 
the lepton-tagged event sample yield consistent results. 

The background-subtracted, efficiency-corrected Stt mass spectrum from the all-tagged sample 
is shown in Fig. ^(a), with the function corresponding to the nominal fit overlaid. Shown in Fig. ^(b) 



is the ai width F^jjj(s), as defined by Eqn. 11, as well as the contributions from the individual ai 
decay channels considered. The kink associated with the turn-on of oi — > K*K is visible in the 37r 
mass spectrum at ~ 1.375 GeV. 

C. Modifications to the nominal fit function 
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FIG. 9. (a) Background-subtracted, efficiency-corrected Svr mass spectrum from the all-tag 
sample. The solid line shows the result from the nominal fit. The inset shows a magnified view of 
the ai -^ K*K threshold region. In (b), the \/s-dependent ai width Tll^[s) is plotted as a function 



of s, along with the separate contributions from the channels considered, as given in Eqn. 11. The 
labels T'^^^^'^ and ri^"*'''^' refer to the contributions from the aj~ -^ ir^ir^ir^ and vr^vr'^vr'^ channels 
respectively. 
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TABLE VI. Results for fits to the Svr resonance shape, with nominal background and efficiency 
correction. The second column gives results based on the nominal fit function. The first error is 
that due to statistics, while the second is the systematic error. The third column gives results from 
the fit including an a'^ contribution, with statistical errors only. The derived quantity B{K*K), the 
ai — > K*K branching fraction, is also shown assuming this is the only amplitude accounted for by 
the threshold function labeled Tk*k in Fig- |9|(b) (see text). The quantity cp^ denotes the phase of 
the a'l amplitude relative to that of the ai . 

Fit Parameter Nominal Fit Fit with a'^ 

moai (GeV) 1.331 ± 0.010 ± 0.003 1.330 ±0.011 

Toai (GeV) 0.814 ±0.036 ±0.013 0.814 ± 0.038 

'ja, (K*K) 3.32 ± 0.26 ± 0.04 3.72 ± 0.45 

B{K*K) {%) 3.3 ±0.5 ±0.1 4.0 ±1.0 

\k\ 0.053 ±0.019 

(f)^/TT — 0.10 ±0.22 

X^/ndof 39.3/41 28.9/39 



1. Finite size of the oi meson 

Like the Dalitz plot distributions, the Sir mass spectrum contains some sensitivity to the 
parametrization of the ai as being point-like or of finite size. Here the sensitivity depends also 
on the treatment of the -y/s-dependence of ai mass. 

As in the substructure analysis, we re-fit the data multiple times, stepping through a range of 
values for the ai meson size parameter R. The results from these fits are given in Appendix |y. For 
both the all-tag and lepton-only-tag samples, the best fits were obtained with values of R between 
1.2 and 1.4 GeV~^, depending on whether the oi mass was treated as a running or constant mass. 
This agrees well with the substructure analysis which favors a value of R of 1.4 GeV~^. This value 
of R corresponds to an r.m.s. radius of the oi of 0.7 fm. It is interesting to note that this is similar 
to the value employed by Isgur, Morningstar and Reader citelMR in their analysis of the 37r line 
shapes from the DELCO |2|, MARK II ^ and ARGUS [||] experiments. However, the statistics 
of the present sample are not sufficient to determine the necessity of including the Fji{k) form 
factor in the parametrization of the hadronic current. 

Despite this, uncertainty on this issue represents a substantial source of model dependence with 
regard to the ai resonance parameters. Tables VIII and IX in Appendix y demonstrate this point. 



As an example, for R = 1.2 GeV with a constant oi mass, we find niQa^ = 1.285 ± 0.007 GeV 
and Foai = 0.619 ± 0.021 GeV (statistical errors only). 

2. Running of the ai mass 

As indicated above, we have also performed fits with the y^-dependence of the oi mass com- 
puted according to Eqn. g. This can only be done in models with non-zero values for the ai size 
parameter R. The results, also given in Appendix ^, indicate that slightly better fits can be ob- 
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tained using a running ai mass. However, since satisfactory fits are obtained with a constant mass, 
we conclude that the present data sample is not sensitive to the running of the ai mass. 

3. Inclusion of an a'i{1700) admixture 

Despite the goodness of the fit to the nominal model, the data above 1.575 GeV show an excess 
relative to the fit function in Fig. 0. This region is where contributions from interference with an 
a'l meson with mass around 1.7 GeV might appear. We have performed various fits allowing k in 
Eqn. to float. The results from one such fit are given in the last column of Table Vl and plotted 



in Fig. |^(a). In this fit, we have used m,oa' = 1.700 GeV and Tq^' = 0.300 GeV. We have also fixed 
the a'i{p'7r) coupling to be equal to the a'l^pir) coupling to determine F^' (s) as shown in Fig. [lO|(b). 
This is an ad hoc choice, however the fit is relatively insensitive to the parametrization of F„/ (s). 

The x^ ^OT this fit is 28.9 for 39 degrees of freedom. The a'^ contribution has a significance 
of 2.8(T, with \k,\ = 0.053 it 0.019 and a phase </>« consistent with zero. If we take (J)k to be zero, 
the resulting fit yields |«;| = 0.060 it 0.017. To test whether backgrounds in the all-tag sample 
are infiuencing this result, we have also fit the spectrum from just the lepton-tag events (with (pf^ 
floating). This fit yields consistent results, with \k\ = 0.05 it 0.04. 



The analysis of DELPHI |12| yielded a large amplitude for the a'l contribution, with |k| in the 
range of 0.50±0.06 to 0.75±0.06, depending on the model used. The corresponding values obtained 
here are smaller by an order of magnitude, and are much less significant statistically. However, the 
models used by DELPHI include neither the contributions to F"oj(s) from the opening of the K*K 
channel, nor those associated with the isoscalar channels. In particular, the K*K amplitude has a 
significant effect at large values of -y/i where the spectrum is most sensitive to the presence of an 
a']^. Omitting the K*K amplitude in our fit, we obtain a value of |fi;| approximately twice as large 
as our nominal value. Finally, the DELPHI anaylsis involves a simultaneous fit to the 37r mass 
spectrum and the Dalitz plot projections at large -y/s, employing assumptions for the substructure 
composition that differ from what we have determined (for the ai) in our substructure analysis. 
In summary, direct comparison of our value of n with that from the DELPHI analysis is not 
meaningful. 

4- Opening o//o(980)7r threshold 

Since the impact of the opening of K*K threshold appeared to be significant, we also performed 
fits including the opening of /o(980)7r threshold. This would include contributions through both 
/o -^ vrvr and KK channels. The /o — > vr'^vr'^ channel would be expected to contribute to the 
substructure in our -k^tt^tt^ sample, but the other modes would not do so. 

Marginal improvements in the fit quality were obtained only for those models in which the oi 
mass was run according to Eqn. |^. In these cases the /o(980)7r contributions to the total width 
of the oi were present at the 2a level or less. Fits including this contribution are discussed in 
Appendix ^ 

To summarize, we find no evidence for contribution from the opening of the oi -^ /o(980)7r 
channel in either the substructure fits or the Stt mass fits. However, other scalar 1 = mesons [a 
and /o(1370)] are needed to provide a good description of the substructure. This observation may 
have some bearing on the interpretation of the /o(980) as something other than a qq meson, as 
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FIG. 10. (a) Background-subtracted, efficiency-corrected Svr mass spectrum from tlie all-tag 
sample. The solid line shows the fit result for the model in which an a[ admixture is allowed to 
contribute. In (b) the assumed pn and p'tt contributions to the a'^ width in this model are shown, 
along with the contributions to the oi width from Fig. pi for comparison. 



has been frequently speculated (see for example Ref. |Q]). Given the theoretical and experimental 
complexity, we cannot comment on this issue except to note that the non-observation of the /o(980) 
in ai decay is not inconsistent with an exotic interpretation for this state. 
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D. Systematic errors 



The sources of the systematic errors shown in Table VI are just those associated with back- 
ground subtraction and acceptance. The background errors are estimated by repeating the fit 
after separately varying the amount of each of the backgrounds being subtracted. We vary the 
background fractions of modes with two real vr'^'s by three times the uncertainty on their branch- 
ing fractions. In the case of the fake-Tr*^ background, the subtraction is varied by three times the 
statistical error of the vr^vr*^ side band sample. Finally, the \/i-dependence of the acceptance as 
determined from Monte Carlo events is parametrized; these parameters are then varied by three 
standard deviations to estimate the associated uncertainty. 



The systematic errors shown in Table VI are dominated by the errors due to the background 
subtraction. As a check, we have performed fits allowing the separate background normalization and 
acceptance correction functions to float, subject to constraints added to the x'^ on the magnitude 
of their deviations from nominal. The changes in fit parameters observed in these fits are small 
relative to the quoted errors. The deviations of the correction functions from nominal are also small 
in this fit. As in the substructure fits, we do not assign a systematic error for model dependence. 



E. Discussion 

As can be seen from Fig. |9|(a), and from the x^ for the nominal fit shown, the model described 
in Sec. [IIIB , with the assumptions listed in Sec. [VIA , provides a good description of the data. The 



oi pole mass and width are determined precisely, however their values depend significantly on the 
model and input assumptions. 

The results demonstrate the importance of including the opening of K*K threshold. Fits per- 
formed without this contribution to Tl^^{s) yield large x^ values. The obtained ai -^ K*K branch- 
ing fraction (3.3 it 0.5)% corresponds to a r branching fraction of B{t~ -^ Ut-o-i -^ Vr[K* K]~) = 
(0.16±0.03)%. For comparison, multiplying the directly measured r^ -^ VrK*^K~ branching frac- 
tion § by a factor of two (to account for K*- K^) gives B{t- -^ Ur[K*K]-) = (0.42±0.08)%. The 
apparent shortfall in our measurement is not surprising since the i/t-[K*K]^ final state is expected 
to receive contributions from both vector and axial vector hadronic currents. 

In the high-mass region above 1.575 GeV, the data appear to be systematically high relative 
to the fit function. To understand whether this is an experimental effect, we have modified the 
background and efficiency corrections within reasonable limits. As described in the previous section, 
we have also performed fits in which background and efficiency correction functions are allowed to 
float. Neither of these approaches significantly improves the fit in this region. Fits with non-zero 
values of Ri, with and without a running oi mass, and/or including the /o(980)7r threshold also 
fail in this respect. 

Including a contribution from an a'l meson, however, visibly infiuences the shape of the 3tt mass 
spectrum in the high-mass region, as shown in Fig. [in(a). As noted earlier, the presence of an a'^ 
may also be consistent with an enhanced d-wave contribution from yo(1450)7r relative to s-wave, as 
observed in the substructure analysis. We have not evaluated systematic errors associated with the 
a[ contribution determined from the 37r mass spectrum fits, since these are likely dominated by 
uncertainties associated with the modeling of the oi line shape. We have performed other fits with 
an a'l, sampling the range of model variations described above. The statistical significance of the 
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a'l contribution is typically 2-3 o" in these fits. We conclude that more data is needed to establish 
whether the a'l is present. 

VII. SUMMARY AND CONCLUSIONS 

To summarize, we have presented a detailed model-dependent analysis of hadronic structure 
in the decay t~ -^ Vt-tt^tt^tt^ using data obtained with the CLEG II detector. This decay mode 
represents a unique source of information on the 7 = 1 axial vector meson sector, an area of hadron 
spectroscopy which is difficult to access cleanly via other production mechanisms. In this analysis 
we have derived successful descriptions of both the ai line shape and the substructure present in 
its decay to three pions. The most significant result is the observation of large contributions to 
the substructure from intermediate states involving the isoscalar mesons a, /o(1370), and /2(1270). 
With this, our data also provides new input on the complicated / = scalar meson sector: we 
observe some sensitivity to the properties of the a meson, for example. More generally, significant 
progress towards a satisfactory description of Dalitz plot distributions in t~ -^ i^^fS-Tr]" decay has 
been achieved. This is supported by the observation that our characterization of substructure in 
the i'-riT~TT^TT^ uiodc also provides a good description of substructure in the v-rir^ tt^ tt^ data. 

Using the results from the substructure fits to infer the \/s-dependence of the oi width, we 
have determined the ai meson resonance parameters. We obtain m-oai = 1.331 it 0.010 it 0.003 GeV 
and Toai = 0.814 it 0.036 ± 0.013 GeV, although these values depend significantly on the details 
of the model used to fit the Sir mass spectrum. For example, taking the meson size parameter 
to be i? = 1.2GeV~ instead of zero, we find rriQa^ = 1.285 it 0.007 (stat.) GeV and Foai = 
0.619 lb 0.021 (stat.) GeV. Such model dependences are not reflected in the quoted systematic 
errors. We also find a significant contribution to the y^-dependence of the ai width associated 
with the opening of the ai — > K*K decay channel at high values of ^/s. 

We have investigated the possibility of an additional contribution to the 37r mass spectrum 



from a radially excited a'l meson, as suggested by the analysis of DELPHI |11,12|. This is also 
suggested by our data, which show an apparent excess of events at large 2>tt mass relative to various 
fits without an o'^ component. The data are better described with an a'l contribution, though at 
a level below that reported by DELPHI. The model used in our analysis differs substantially from 
that analysis, with regard to treatment of the K*K threshold and the substructure in the ai — > 37r 
channel. We have not assessed the impact on v-j- mass studies of effects associated with the complex 
substructure in ai -^ Stt decay or the apparent distortions in the Sir mass spectrum caused by the 
K*K and possible a'^ contributions. However, careful consideration of such effects in the course of 
these analyses should improve the reliability of the ensuing tau neutrino mass constraints. 

We have also obtained a precise determination of the signed r neutrino helicity h^^ = 
—1.02 lb 0.13 (stat.) lb 0.03 (syst. + model) when this quantity is left as a free parameter in the 
substructure fits. As has been noted in earlier measurements of this quantity Q, accurate pa- 
rameterization of the substructure is important for obtaining an unbiased measurement. With the 
improved understanding of this substructure, this result provides unambiguous evidence for the 
left-handedness of the r neutrino. 

We have addressed several additional issues pertaining to the characterization of axial vector 
meson decay dynamics. For example, the data show limited sensitivity to the finite size of the ai 
meson, both in the substructure and the 37r mass spectrum fits. With the parametrization of the 
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associated form factor used here, we find that both analyses favor an r.m.s. ai radius of around 
0.7 fm. As with the question of a possible a'l contribution, a definitive conclusion on this issue 
requires additional data. We have looked for indications of non-axial- vector contributions to the 
Stt substructure via the 7r'(1300) resonance, and have placed upper limits on the r decay rate to 
this state. Detailed analyses of the higher-statistics t~ — > ^'T-vr^vr+vr^ data should shed additional 
light on these and other issues. 

Although the quantitative results presented in this article are strongly model-dependent, they 
describe successfully the qualitative features of the data. However, further insight can be gained 
from a quantitative model-independent analysis of the data, such as that proposed by Kiihn and 
Mirkes [29|. We are presently pursuing such an analysis and will report the results in a separate 



article. 
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APPENDIX A: PARAMETERIZATION OF SUBSTRUCTURE IN r" -^ i^^n'Tv^Tv^ 

To parameterize the model used to fit for tlie substructure in r~ — > i/T-vr^Tr^vr^, we follow many 



of the conventions used in the KORALB/TAUOLA [17| tt Monte Carlo generator. We first denote the 
four-momenta of the pions by 

pi = four-momentum of tt^, 

P2 = four-momentum of 7r2, and 

P3 = four-momentum of vr^. (Al) 

We then define the quantities qi = P2—P3, Q2 = Ps—Pi, Q3 =Pi—P2, and a = pi+P2+P3- Ignoring 
for now the resonant structure of the Svr system, the general form for Svr hadronic current J^ as 
defined by Eqs. Q and ^ can be written as: 

J" = T"-" [cm^Fi + C2q2uF2 + cm^Fs] + c^a^^Fi + c^ie^'^P'^pi^pipP^aF^ , (A2) 

where the q are complex scalar coefficients, and T'^'^ = g'^'^ — a'^a'^ /a?. 

The form factors Fi contain the description of the low energy QCD phenomena we are studying. 
Owing to Lorentz invariance, they depend only on Lorentz scalars. The terms containing i<i, F2 
and i^3 are associated with the axial vector contribution. The term containing F^, is not needed 
since it can be absorbed into the terms containing Fi and F2. We write it explicitly here to 
make connection with form factors associated with specific resonant substructure. F4 is the scalar 
form factor, and F^ is the G-parity violating vector form factor. Neither of these are expected to 
contribute significantly in t~ — > ^'T-[37r]~, hence we generally set these terms to zero, except where 
noted. 



1. Model for the form factors 



The form factors as defined in Eqn. A2 do not have a simple correspondance with those that can 
be associated with specific resonant contributions to the hadronic current. Here we give the ansatz 
for the amplitudes j^ (as defined in Eqn. |3|) of the hadronic current in the decay r~ -^ Vt-tt^'k^tt^ 
used in the substructure fits: 
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jf^ —'j^fiiy 



qi,BP{si)FR,{ki) - q2uBP{s2)FR,{k2) 



j^^T^u 



qi,BPisi)FR,ih) - q2uBP,{s2)FR,{k2) 



j^^rj.^U 



Qi, {aqi) BP{si)FR,{ki) - Q2. (0^2) 5f (82)^^3 (/^s) 



j^=T^'^ 



h. (aqi) BP{si)FR,{ki) - Q2. (092) BP{s2)FR,{k2) 



(A3) 



J^=T^'' 



qsu (aqs 



o-v — hsu 



(hsa) 



S3 



iQsqs] 



B%{ss)FR,ik3 



J^=T^'' 



Q3uB^{s3)FR^{k3 



ji^^J.^u 



Q3uBUss)FR,{k3 



The kinematic factors appearing in above are defined as follows: 

hi = Pj + Pk with {i^ j ^ k ^i) 

n; J — '^i Pi : 

■ — h^ — 2 

and the decay momenta are given by 



k 



Ji- 


- {y/sl + mi)'^j [s- 


- {y/s^i-mifj 


I 


2^s 




._h 


- {rrij + mfc)^j [si 


- {rrij - rukfj 



Z-^Si 



with {i i^ j i^ k ^ i) 



(A4) 



(A5) 



(A6) 



For the complex couplings /?j, we specify /?i = 1, and thus we determine the remaining couplings 
relative to the first amplitude (pvr, s-wave). The Breit-Wigner functions for the intermediate states 
are 



B^{si) 



m: 



OY 



{mlY - Si) - imor'i^^'^isi) 



^^'^-)-^^(|)""T 



(A7) 
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The parameters moY, Tg , k'^ and /cq are the nominal mass, nominal width, the decay momentum, 
and the decay momentum at Si = m-gy, respectively. See Table | for a summary of the resonance 
parameters of the intermediate states as used in the substructure fits. The ansatz used for the form 
factors Fr^ is 



FnAh 



exp( 



I^m- 



(A8) 



2. Connection with reduced form factors 



With the ansatz J'^ = J2i=i PiJi^ the reduced form factors Fj of Eqn. A2 are given by: 
Fi= (3iB^{si)FR,{ki) 



+ p2BJ^,{si)FR,{ki 



ms3 - ml) - (si - ml))BJ^is2)FR,{k2) 
iMis3 - ml) - (si - ml))Bp{s2)FR,{k2) 



+y. ^"' - "^^ + ^^^?"' + ^"^ - ''^ B-is,)FR,ik,) 



6S3 



+^f3eB^{s3)FR,{k3 



+ PjB%{ss)FR,{k3 



(A9) 



F^ 



^^3 i{s2 - ml) - (S3 - ml))B';{si)FR,{k^) + ((53 - ml) - (si - ml))B'; {s2)FR,{k2) 



+ 3/34 [((S2 - ml) - (S3 - ml))Bp{si)FR,{ki) + ((S3 - ml) - (si - ml))BP,{s2)FR,{k2) 



-2/35((5i - ml) - (S2 - ml))Bl{s^)FR,{k,) 







(AlO) 
(All) 



F2 is given by Fi under an interchange of indices 1 <-> 2, with a relative minus sign between Fi and 
F2 required (i.e., the coefficients Cj Eqn. |A^ follow C2 = — ci). 



APPENDIX B: SUBSTRUCTURE FITS EXCLUDING ISOSCALARS 



To assess the statistical significances of individual amplitudes in the nominal substructure fit 



summarized in Table III, we successively repeated the fit, each time with one amplitude omitted. In 
view of the large contribution from amplitudes involving isoscalar mesons, we have also performed 
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a fit in which all three of these have been omitted. The Dalitz plot projections from this fit are 



plotted in slices of ^/s in Figs. 11 and |12|, overlaid on the data distributions. The agreement with 
the data is visibly worse than that seen in the nominal fit with all amplitudes (see Figs. |3| and Q). 
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FIG. 11. Dalitz plot projections: distributions in squared 7r~-K^ mass si and S2 (two entries per 
event), plotted for slices of ^/s as in Fig. ^. The data are represented by the filled points. The solid 
line is the result from a modified version of the nominal fit in which amplitudes involving isoscalar 
mesons (crvr, /qtt and /27r) have been omitted. 



APPENDIX C: RESULTS FROM FITS WITH FINITE MESON RADII 

Treatment of the ai as being point-like results in unphysical behavior of the running width 
r"oi(s) at large values of s. The effect of the form factor FR^{ki) defined in Eqn. ^ is to damp out 
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FIG. 12. Dalitz plot projections: distributions in squared vr'^vr'^ mass S3 (one entry per event), 
plotted for slices of ^/s as in Fig. |^. The data are represented by the filled points. The solid line is 
the result from a modified version of the nominal fit in which amplitudes involving isoscalar mesons 
((Tvr, /oTT and /27r) have been omitted. 

this behavior, although the Gaussian form itself is somewhat ad hoc. This form factor affects both 
the substructure analysis and the Svr mass spectrum analysis. Although the data at present do not 
show much sensitivity to the value of the meson size R, we note that both analyses prefer values 
of i? ~ 1.4 GeV^ . In this appendix we report the results from fits using non-zero values for R. 
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1. Substructure fits 



In Table |VI]| , we give the results from the substructure fit with Ri 
are in qualitative agreement with those from the nominal fit. 



1.4 GeV . These results 



TABLE VII. Results from substructure fit with meson radii Ri = 1.4 GeV ^. 







Signif. 


m 


•/•ft/TT 


B fraction (%) 


p 


s-wave 


- 


1.00 


0.0 


69.11 


/9(1450) 


s-wave 


0.6a 


0.03 ± 0.06 


0.92 ±0.58 


0.02 ± 0.08 


P 


d-wave 


5.4a 


0.36 ± 0.09 


-0.12 ±0.11 


0.28 ±0.13 


p(1450) 


d-wave 


3.7a 


0.78 ±0.28 


0.41 ±0.14 


0.28 ±0.20 


/2(1270) 


p-wave 


4.6a 


0.66 ±0.14 


0.64 ±0.10 


0.11 ±0.05 


a 


p-wave 


7.9a 


2.13 ±0.27 


0.27 ±0.03 


14.31 ±3.66 


/o(1370) 


p-wave 


6.1a 


0.77 ±0.13 


-0.47 ±0.05 


3.69 ±1.21 



2. Fits to Mstt 



Fits to the Svr mass spectrum were performed with non-zero values for Ri with the assumptions 
of constant and running ai masses. These studies differed from the nominal fit in Fig. M in two 
additional ways. First, the fitting range included two additional bins in the high mass region, 
extending up to 1.775 MeV. Second, so as to account for possible systematic effects, the acceptance 
and background corrections were allowed to vary within reasonable amounts by adding terms to 
the x^ constraining their deviations from the nominal assuming these to be Gaussian distributed 
(see Section VI D| ). 



Results for fits assuming constant and running oi masses are presented in Tables VIII and [X 
respectively. Similar trends are observed in fits to just the lepton-tagged sample, as well as in fits 
that include an a'^ contribution. The final column in the tables, labeled y/x^, denotes the square 
root of the contribution of the acceptance and background correction constraints to the total x^ of 
the fit. In both constant and running mass scenarios, the oi pole width is strongly affected by the 
presence of the form factor accounting for the size of the oi meson, and by the value of R. 

The preferred fits, with R = 1.2 GeV~ for the constant mass case and R = 1.4 GeV^ for the 
running mass case are plotted in Figures |l^( a) and (c), respectively. We have also performed these 
fits including turn-on of the /o(980)7r channel, the results of which are plotted in Figures 13(b) 
and (d) for the two cases. 

The parametrizations of the -y/i-dependence of the ai mass entering the fits shown in Figs. p!^(c) 
and (d) are plotted in Fig. |l^. The overall mass shift function depends on the relative amplitudes 
for the K*K and (in Fig. |T3|(d)) /o(980)7r channels which are fit parameters, however the shapes of 



the contributions to mai{s) from these channels are determined as described in Section IIIB . The 



effect of successively including thresholds is illustrated by the dotted curve in Fig. 14(a), in which 
an ad hoc /o(980)7r contribution is added, assuming unchanged couplings to the other channels. As 
noted earlier, the effect is to flatten the dependence of mai{s). 
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TABLE VIII. Results of Sir mass spectrum fits with different input values of oi size parameter 
R, assuming constant ai mass. 

R X^/ndof rriQa^ Toai B{K*K) ^^ 

[GeV-^] [GeV] [GeV] [%] 



0.0 


40.9/43 


1.330 ±0.011 


0.811 ±0.042 


3.3 ±0.6 


0.32 


1.0 


38.9/43 


1.289 ± 0.008 


0.653 ± 0.025 


3.4 ±0.6 


0.32 


1.2 


38.6/43 


1.285 ± 0.007 


0.619 ±0.021 


3.5 ±0.6 


0.31 


1.3 


39.3/43 


1.279 ± 0.006 


0.597 ±0.019 


3.4 ±0.7 


0.33 


1.4 


39.9/43 


1.274 ± 0.006 


0.578 ±0.017 


3.2 ±0.7 


0.37 


1.5 


42.9/43 


1.269 ±0.005 


0.558 ±0.015 


3.0 ±0.7 


0.38 


1.6 


45.5/43 


1.263 ±0.005 


0.538 ± 0.014 


2.8 ±0.7 


0.44 



TABLE IX. Results of 37r mass spectrum fits with different input values of ai size parameter 
i?, assuming running ai mass. 

R X^/ndof rnoaj Loai B{K*K) Vxl 

[GeV-^] [GeV] [GeV] [%] 



1.2 


39.6/43 


1.298 ± 0.007 


1.200 ±0.100 


3.3 ±0.7 


0.77 


1.3 


39.3/43 


1.294 ± 0.006 


0.822 ± 0.047 


3.3 ±0.7 


0.78 


1.4 


36.7/43 


1.288 ± 0.006 


0.667 ±0.031 


3.3 ±0.7 


0.71 


1.5 


41.9/43 


1.284 ± 0.006 


0.558 ± 0.021 


3.2 ±0.7 


0.79 


1.6 


54.4/43 


1.279 ± 0.005 


0.493 ± 0.016 


3.0 ±0.7 


0.77 
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FIG. 13. Background-subtracted efficiency-corrected Svr mass spectrum (points with error bars) 
witli fit functions overlaid. The solid curves in plots (a) and (b) represent the fit function with 
R = 1.2 GeV~ , assuming a constant mass, excluding (a) and including (b) the effect of the 
/o(980)7r threshold on Tll^{s). The solid curves in plots (c) and (d) represent the fit function with 
R = 1.4 GeV~^, assuming a running mass. The effect of the /o(980)7r threshold on Tll^{s) and 
mai(s) is excluded in (c), and included in (d). The dotted curves represent the corresponding fit 
functions without the deviations in the background subtraction and efficiency correction returned 
by the fits. 
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FIG. 14. The solid curve in plot (a) shows the running mass used in the fit shown in Fig. ^(c). 
The dotted curve illustrates the effect of adding the /o(980)7r threshold, with r37r(s) and Tk*k 



unchanged. Plot (b) shows the running mass used in the fit shown in Fig. 13(d), where the relative 
K*K and /o(980)7r amplitudes are determined from the fit. 



Satisfactory fits are obtained with both constant and running ai masses. Functions assuming 
a running mass yield fits with smaller x^ values than those with a constant mass. However, the 
running mass fits also prefer a larger distortion of the background and acceptance corrections, as 
indicated by the larger values of 



Xc in Table IX and by the deviations of the dotted curves in 



Figs. |T^(c) and (d) from the corresponding solid curves. 
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